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Future multifrequency microwave background experiments with arcminute resolution and micro- 
Kelvin temperature sensitivity will be able to detect the kinetic Sunyaev-Zeldovich effect, providing 
a way to measure radial peculiar velocities of massive galajcy clusters. We show that cluster peculiar 
velocities have the potential to constrain several dark energy parameters. We compare three ve- 
locity statistics (the distribution of radial velocities, the mean pairwise streaming velocity, and the 
velocity correlation function) and analyze the relative merits of these statistics in constraining dark 
energy parameters. Of the three statistics, mean pairwise streaming velocity provides constraints 
that are least sensitive to velocity errors: the constraints on parameters degrades only by a factor 
00 ' of two when the random error is increased from 100 to 500 km/s. We also compare cluster veloc- 

L~:> ' ities with other dark energy probes proposed in the Dark Energy Task Force report. For cluster 

I velocity measurements with realistic priors, the eventual constraints on the dark energy density, 

■ the dark energy equation of state and its evolution are comparable to constraints from supernovae 

5^ ' measurements, and better than cluster counts and baryon acoustic oscillations; adding velocity to 

other dark energy probes improves constraints on the figure of merit by more than a factor of two. 
<^ . For upcoming Sunyaev-Zeldovich galaxy cluster surveys, even velocity measurements with errors as 

' large as 1000 km/s will substantially improve the cosmological constraints compared to using the 

'/^ ' cluster number density alone. 

<N : 

. PACS numbers: 95.36.-l-x, 98.65.Cw, 98.80.-k 

6 ■ 

^ ■ I. INTRODUCTION 

. . . . 

Historically, the primary goal of cosmology has been determination of cosmological parameters describing the 
overall properties of the universe. This quest has advanced greatly in the past decade with the precise measurement 
' of microwave background temperature fluctuations putting sharp constraints on many parameters [l|, Pushing 
J|T 1 cosmological parameter determination to ever-greater precision might have been an academic pursuit, except for the 
' surprising discoveryof the universe's accelerating expansion [sj, 01 , coupled with the discrepancy between the geometry 
of the universe [1, Q and its mass density 0, Ullj ■ The likely implications of the inferred "dark energy" 

for fundamental physics are fueling a wide array of next-generation cosmological experiments. Current measurements 
I can constrain the dark energy density and its current equation of state at around the 10% level from a combination of 
I • WMAP 3-year data and the Supernova Legacy Survey 1-year data [isj , for example. However, current experiments 
J — I can only constrain dark energy density and its equation of state at the present epoch. We also need to quantify the 
. evolution of dark energy with redshift, for this determines whether the dark energy is fundamentally a dynamic 
' (evolving scalar field) or static (cosmological constant) entity. In recent years, a quartet of methods has emerged as 
• i-H , the most discussed for constraining dark energy redshift evolution: weak lensing by large-scale structure, primordial 
' baryon acoustic oscillations (BAO) observed as a feature in the matter power spectrum at low redshifts, the distance- 
JH [ redshift relation measured via SNIa standard candles, and the redshift evolution of galaxy cluster counts detected via 
■ - - i the Sunyaev-Zeldovich Effect. The relative merits of these probes were considered in detail by the recent report from 
the Dark Energy Task Force 14]. Each of these probes suffers from different sources of systematic error: cluster counts 
are subject to uncertainty in the mass-SZ relation[l^; weak lensing and BAO suffer from uncertainty in modelling 
baryon physics and nonlinear effects [l^, [3 i ^^'^ SNIa distance mesurements require knowledge about the extent 
to which the supernovas serve as standardizable candles over a range of redshifts fi3\ . It is clearly important to probe 
cosmology through multiple techniques to check consistency between individual probes. 

This paper addresses another approach to constraining dark energy which so far has received comparatively little 
attention: the line-of-sight peculiar velocities of galaxy clusters. A moving galaxy cluster will induce a nearly- 
blackbody shift in the distribution of the microwave photons passing through the cluster, due to Compton upscattering 
of the photons by hot electrons in the cluster gas (20| . The temperature shift, known as the kinematic or kinetic 
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Sunyaev-Zeldovich (kSZ) effect, is proportional to the line-of-sight momentum of the cluster gas (being linearly 
proportional to both the optical depth for Compton scattering and to the line-of-sight velocity of the cluster with 
respect to the microwave background rest frame), while being independent of the cluster gas temperature. It is 
substantially smaller than the more familiar thermal Sunyaev-Zeldovich effect [2l|: a typical large cluster with a 
thermal SZ distortion of 100 /iK at the frequency with the largest distortion may have a kSZ signal of 5 to 10 /iK for 
typical cluster velocities. The kSZ effect can be thought of as essentially a Doppler shift due to cluster motion. Like 
the related thermal SZ effect, the kSZ effect has the remarkable property that its imprint in the microwave background 
which we observe today is independent of the cluster's redshift, making it potentially an excellent probe of cosmology. 
If we can reliably measure the kSZ effect in galaxy clusters, we expect the line-of-sight velocity error for an individual 
cluster to be largely independent of the cluster redshift, in marked contrast to galaxy peculiar velocity surveys. 

The cosmological velocity field on cluster scales, arising solely from the effects of gravitational instability in the 
universe, is a potent probe of structure formation. It is thus also potentially a strong probe of dark energy, whose 
properties affect the rate of structure growth. Here we study the feasibility of probing dark energy parameters using 
peculiar velocity of galaxy clusters obtainable through detection of kinematic Sunyaev-Zeldovich effect. The small 
amplitude of kSZ distortions, combined with the need to separate this small blackbody signal from several other 
larger signals with various spectra, makes detection challenging. Several studies have shown that in the absense of 
any foreground contamination, it would be possible to measure cluster velocities with reasonable accuracy 100 
to 150 km/s) through multi-frequency SZ measurements with arcminute resolution, combined with X-ray foUowup 
[2^ [23j . Even in the presence of point source contamination, it is possible to measure velocities with an accuracy 
of perhaps 200 km/s through multifrequency measurements with arcminute resolution [24| and sufficient integration 
time. Internal motions of the intracluster medium give an irreducible random error of around 100 km/s (25j . 

Upcoming SZ measurements like ACT [2^,[23 and SPT (2^ are designed to detect large numbers of clusters through 
their SZ signatures. The ACT collaboration foresees maps of sufficient raw sensitivity to measure the kSZ effect in 
many clusters, making detailed studies of the cosmological impact of future kSZ measurements timely. Some recent 
work has shown that the kSZ correlation function will put significant constraints on the dark energy equation of 
state [2^, and cross-correlation of the kSZ signal with the galaxy density can constrain the redshift evolution of the 
equation of state [13] ■ Cluster velocities alone can be used to constrain the matter density of the universe [3l|, [13] , 
the primordial power spectrum normalization [32| , and the dark energy equation of state [32j ] . 

The goal of this paper is twofold. Following up on our initial study [321 ] . we study the accuracy of theoretical models 
of cluster velocity statistics by comparing with numerical simulations, address error analysis in greater detail, and 
compare the relative merits of various velocity statistics in constraining dark energy parameters; we also compare 
cluster velocities with the Dark Energy Task Force methods. We find that all three velocity statistics considered 
here can be computed using the halo model within likely measurement uncertainties. We then use a Fisher matrix 
calculation to compare the power of various velocity statistics as dark energy probes over a range of velocity errors. 
Remarkably, for a sufficiently large velocity catalog, the dark energy parameter constraints degrade only by a factor of 
two when the velocity errors increase by a factor of five. Comparing with other dark energy probes, cluster velocities 
from a large survey can provide dark energy constraints that are comparable to weak lensing and supernovae and 
a factor of two to three better than cluster counts and BAO. Combining cluster velocities with other dark energy 
probes improves the total constraint on the dark energy density by 10-15% and the Dark Energy Task Force Figure 
of Merit by a factor of 1.4 to 2.5. Cluster velocities can be competitive with other proposed techniques for probing 
dark energy, with completely different systematic errors. 

Throughout this paper, we assume a cluster velocity catalog with some normal velocity error; we consider errors 
from 200 km/sec to 1000 km/sec, representing a range from optimistic to conservative based on current experimental 
sensitivities and anticipated astrophysical complications. Using this range of errors, we then evaluate the statistical 
constraints on dark energy parameters, assuming a cluster catalog with a given number of cluster velocities. In 
practice, constraints from cluster velocities may well be dominated by systematic, rather than statistical, errors, like 
all other methods of probing dark energy. When analyzing real data to constrain dark energy, understanding these 
systematic errors is obviously crucial to getting the right answers. For our purposes here, we only aim to evaluate the 
statistical power of cluster velocities to constrain dark energy; we thus ignore systematic errors, keeping in mind that 
any results here require an understanding of all relevant systematics to be realized in practice. Note that systematic 
errors will tend to bias parameter constraints but will not generally change the size of the statistical errors significantly. 
A discussion of various relevant systematics is given in the last section of the paper; we will address this issue in more 
detail elsewhere. 

This paper is organized as follows. Section II gives theoretical approximations of various velocity statistics computed 
using the halo model; Section III studies the accuracy of the halo model expressions by comparing with simulations. 
Section IV discusses various sources of errors for each of the statistics and presents analytic expressions for the 
errors; detailed derivations of these expressions are given in three Appendices. Using these expressions for the values 
of the velocity statistics and their errors in hypothetical surveys of given sky area and velocity errors. Section V 
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uses standard Fisher matrix techniques to compute constraints on dark energy parameters from the various velocity 
statistics. Section VI then compares the cosmological constraints obtainable from cluster velocities with those from 
the probes analyzed by the Dark Energy Task Force. Finally, Section VII discusses further refinement of the current 
calculations, including correlations between various velocity statistics, extraction of cluster velocities from microwave 
maps, and near-future prospects for kSZ velocity measurements. Throughout, we employ a standard spatially flat 
ACDM model with parameters given by the best-fit WMAP 3- year values as our fiducial cosmology unless otherwise 
noted. 



II. THE HALO MODEL FOR VELOCITY STATISTICS 



To study the potential of galaxy cluster velocity surveys to serve as a dark energy probe, we consider three different 
velocity statistics: the probability distribution function of the line-of-sight component of peculiar velocities n„; the 
mean pairwise streaming velocity Vij{r), which is the relative velocity along the line of separation of cluster pairs 
averaged over all pairs at fixed separation r; and the two- point velocity correlation function {viVj){r) as a function 
of separation r. In the halo model picture of the dark matter distribution (33l. [3^. these quantities can be written 
as the sum of the contribution from one-halo and two-halo terms. However, we are interested only in very massive 
clusters {M > IO^^Mq), so the one-halo term can be neglected. 

Here we summarize the halo model ingredients which go into computing the values of these velocity statistics for 
given cosmological models. Define moments of the initial mass distribution with power spectrum P{k) by [ssj 

1 f°° 

cr^(m) = — dk(k^+^^)P(k)W^(kR(m)) (1) 

when smoothed on the scale R(m) = {Sm/AnpoY^"^ with the top-hat filter W{x) = 3[sin(2;) — a;cos(a;)]/a;'^, and po 
the present mean matter density. The spherical top-hat halo profile is adopted for simplicity. It could be replaced 
by a more realistic NFW profile; however, we are interested in statistics only of the most massive clusters at large 
scales where details of halo profiles make no significant difference. We also write H{a) for the Hubble parameter 
as a function of scale factor, h « 0.7 as the Hubble parameter today in units of 100 km/s Mpc~^, and i?iocai for a 
smoothing scale with which the local background density S is defined. 

The number density of halos of a given mass n{m) is taken as the Jenkins mass function [3^ 

— m,z =0.315-^— -^^exp - 0.61 - In ao m Z?, . 2 
am alnm L J 

This mass function is a fit to numerical simulations of cold dark matter gravitational clustering. The bias factor can 
be written as [131 

= 1 + (3) 

where Da is the linear growth factor at scale factor a, normalized to 1 today, and the critical overdensity (5crit ~ 1.686. 
Since clusters preferentially form at points in space of larger overdensity, the number density of clusters for a given 
mass and formed in a given local overdensity can be written as 38] 

n{m\5) w [1 + h{m)S\ n{m). (4) 

The matter power spectrum P{k) at the present epoch can be well fit through a transfer function as 

P{k) = — (5) 

[1 + [ak + (/3fc)3/2 + {^^kYYY'" 

where a = {QA/T)h-^ Mpc, (3 = {i.Q/T)h-^ Mpc, 7 = {1.7 /T)h-^ Mpc, v = 1.13 and F = [H, The 

normalization B is fixed at large scales by normalizing to the microwave background fluctuation amplitude. 



A. Probability Distribution Function 



The probability p{v \ m, 5, a) that a cluster of mass m located in an overdensity 5 moves with a line-of-sight velocity 
V can be approximated by a normal distribution j38j , 



p{v I m, (5, a) 



27r 



1/2 



1 



ay{'m, a) 



exp 



3v 
(Jy{m, a) 



1 



(6) 
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with the three-dimensional velocity dispersion smoothed over a length scale R{m) given by [4lj 

1/2 

a„(m, a) = [1 + <5(i?iocai)]'^^'''-'^ aH{a)Dj^^^ U - ^74^^^) a) (7) 

and Hi] 

Ai(i?iocai) = 0.6a2(i?io,,i)/'^o(10Mpc/h). (8) 

Following i?iocai is obtained empirically using N-body simulations via the condition cro(^iocai) = 0.5(1 + 

Then the probability density function of the line-of-sight peculiar velocity component at some redshift z is given 
by M 

^ J d'mmn{m\S)p{v\'m,S,a) 
' J dmmn{m\S,a) 

where n(m\S)dm is the number density of halos that have mass between m and m + dm in a region with overdensity 
6. The dependence of these quantities on redshift is left implicit. 

Finally, in order to connect to a readily observable quantity, we write the fraction of clusters that have velocity 
between v and v + 6v as 

nv{v,Sv,a) = / dvf{v,a). (10) 

Jsv 



B. Mean Pairwise Velocity 

The mean relative peculiar velocity Vij (r) between all pairs of halos at comoving separation r and scale factor a can 
be related to the linear two-point correlation function for dark matter using large-scale bias and the pair conservation 
equation [i^ : 



2^. , din Da r^^'°{r,a) 

-—H[a)a , , , — , . 

3 ^ ' dlna H-C^^i°(r,a) 



^,„-(r,a) = H{a)a ^ ^ ,1 \ . (11) 



The two-point correlation function can be computed via 



^haio(^^ a)^^ I dkk sin fcrP(fc)6(l, (12) 



while the two-point correlation function averaged over a sphere of radius r can be written as 

7-j2 rr /•CO 

^-i>-i°(r-,a) = -^ / drr / dkk sin krP{k)bi'^,^ (13) 
r Jq Jq 

where average halo bias factors are given by 



^(g) _ J dmmn{m)b{'m)''W'^[kR{m)] 
J dm mn(m)W-^ [kR{m)] 



"halo — f J_ _„/™Mi/9r;,D^„,M • \^^> 



Direct evaluation of the above expression for mean pairwise peculiar velocity requires knowledge of all three velocity 
components for both halos. In practice, it is only possible to determine the radial velocity components, so we need an 
estimator which depends only on the radial velocities. Consider two clusters at positions and Tj moving with 
velocities and Vj . The radial components of velocities can be written as w[ — Vi ■ Wi and — ij ■ Vj . Following 
p3 |. {vl — Vj) = w^l'f • [fi -I- fj]/2 where r is the unit vector along the line joining the two clusters and f is the unit 
vector in the direction r. Then minimizing gives 

where pij = r ■ (vi + Vj) and the sums are over all pairs of clusters with separation r. 
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C. Velocity Correlation Function 

In addition to the mean relative peculiar velocity between two halos, we can also consider correlations of these 
velocities. Assuming statistical isotropy, the only non-trivial correlations will be of the velocity components along 
the line connecting the clusters and of the velocity components perpendicular to the line connecting the clusters; 
furthermore, these correlations will only depend on the separation r 
radial velocities must be of the form 



only 

m 



\ri — Yj\. Geometrically, the correlation of 



'I'_lcos6'+ ('fii - 



(rf + ) cos 6 — rirj{l + cos^ 6) 
rf + r'j — 2rirj cos 6 



(16) 



where = Vi -Vj is the angle between the two cluster positions; ^±{r) and ^|| (r) denotes the correlations perpendicular 
to the line of separation r and parallel to it, respectively. Including the fact that high-density regions have lower rms 
velocities than randoni patches and allowing the two halos to have different masses, the expressions for correlations 
can be written as [43, WOl 



ao{mi)ao{mj) 2 H{af 
^fx, II (mi, TOj-r) = — ra 



din Da 
din a 



Dl / dkP{k)W[kR{m.,)]W[kR{mj)]K^^\\{kr) (17) 



where 



jijkr) 
kr 



K\\ = joikr) 



kr 



(18) 



with joikr) and ji{kr) the spherical Bessel functions. 

With all the above ingredients, the correlation function for the velocity components perpendicular to the line 
connecting the clusters can be written as 



{vtVj)±ir, a) 



d\nDa 

H(a)a— Da 

a ma 



minimi) f , m-jn(mi) 1 + b(mi)b(mi)f(r) , , , , , 

dm, ' ^ I dm, ' : " h ] W ^1 '^Am,m, r) 19 



where p — J dmmn{m) . Note that the above expression is a slight modification from Eq. (23) of Ref. [i^. The expres- 
sion for the correlation of the parallel velocity components is obtained simply by replacing with ^||. Performing 
the ensemble average yields 



{viVj)j_{r,a) = a^H{a)'^ 



2 ( din Dg 
din a 



D 



1 



1 



/i+e^^'°(r,a)/2 



"l+^halo(^^^)2^2 



where 



(20) 



h = / dkK^{kr)P{k) 



dmmn{m) — \ w[kR{m)\ 
(T_i(m) 



(21) 



dkKj_{kr)P{k) 



dmmn{m)b{m) — \ w[kR(rn)] 
cr_i(m) 



(22) 



Although the above expression holds for both the parallel and perpendicular components, in simulations is 
mostly negative or zero due to the heavy influence of infall at large separations j4j|. However, this anticorrclation is 
not seen in linear perturbation theory or in the halo model, which both predict positive correlation for pair comoving 
separations less than 40 Mpc; for separations larger than 40 Mpc, the theory and simulations are consistent, but 
the parallel component correlation is essentially zero. Given this discrepancy between known analytical models and 
simulations for the parallel correlation function in the region where the signal is non-negligible, we only consider 
(viVj) j_{r, z) in the rest of this paper. 



III. COMPARISON WITH SIMULATIONS 



The statistics computed in the previous section are based on the halo model of structure formation combined with 
linear perturbation theory. Since galaxy clusters are rare objects and their distribution can be described well in 
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FIG. 1: A comparison between the probability distribution function n„ evaluated directly using the Virgo lightcone numerical 
simulation (dotted curve with error bars) and approximated using the analytic halo model formula, Eq. (|10p (solid red curve). 
Error bars are Poisson plus cosmic variance errors for one octant sky coverage. 

the quasi-linear regime of structure formation, we expect that these approximations for velocity statistics should be 
reasonably accurate. Here we verify that they are good approximations to the actual galaxy cluster velocity statistics 
extracted from the the VIRGO dark matter simulation [47|. We use the octant sky survey (PO) lightcone output of 
LCDM cosmology, with erg = 0.9, = 1, Vim = 0.3, JIa = 0.7 and h = 0.7. The maximum redshift of the light cone 
is ^max = 1-46 and the radius of extent is i?max = 3000 Mpc/h. The data is binned in redshift slices of width 6z = 0.2 
from z = to z = 1.4 

The statistics defined in the previous section apply to infinitesimal intervals in redshift. When comparing with 
data binned in redshift, it is necessary to normalize the velocity statistics properly to reflect this binning. We do this 
by averaging the above theoretical expressions for the statistics over a given bin in z to obtain a binned estimator 
of the underlying statistic. Additionally, for the case of the velocity probability distribution function, a realistic 
measurement will provide numbers of clusters in a set of line-of-sight velocity bins. In this case, the relevant statistic 
for comparison becomes the theoretical probability that the velocity of a given cluster is in a particular velocity 
bin; the above expression for probability density in infinitesimal velocity bins must be integrated over the width of 
the velocity bin. This gives the correct relative probabilities between any two velocity bins, but all must then be 
renormalized by a constant factor to enforce the condition that the sum of the probabilities for all bins be unity. 
All comparisons with simulations below use these binned versions of the underlying statistics defined in the previous 
section. 

Figure [T] shows riy in the redshift slice between z = and z = 0.2 both from the simulation and using Eq. (fTU)) 
for the velocity probability distribution function. The analytical model agrees fairly well with the simulation; the 
error bars denote the la errors including Poisson error and errors due to cosmic variance. Error modeling is discussed 
in detail in the next Section. Note that the error bars shown in Figure [T] are for a large future 5000 square degree 
velocity survey (one octant of the sky) . Figures [5] and [3] compare the simulation with Eq. (jlip for the mean pairwise 
streaming velocity (using the estimator Eq. p3|) ) and Eq. (pO)) for the velocity correlation function, respectively. The 
plots shows that the halo model agrees well with the simulated data at separations greater than 30 Mpc/h for velocity 
correlation, and greater than 40 Mpc/h for mean pairwise streaming velocity with a discrepancy somewhat larger 
than 1(7 for r between 30 and 40 Mpc/h. For the velocity probability distribution function, we find a good fit when 
the velocity data is smoothed over a scale of 10 Mpc. The smoothing on this scale reduces the effect of nonlinear 
physics which is difhcult to model semianalytically. 

Figure 3] displays a comparison between the estimated mean pairwise streaming velocity obtained only from 
the radial component of velocity using Eq. (jlSp and the full Vij obtained from all three components of velocity in the 
simulation. For an ideal estimator, these quantities would be exactly the same; the actual estimator in general does 
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FIG. 2: A comparison between the mean pairwise streaming velocity Vij{r) evaluated directly using the Virgo lightcone 
numerical simulation (dashed line with la errors given by the blue dotted lines) and approximated using the analytic halo 
model formula, Eq. ((TTJ (red solid curve). The error range includes Poisson and cosmic variance errors for one octant sky 
coverage, plus random measurement errors of 100 km/s. 
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FIG. 3: Same as Figure [2] except for the velocity correlation {viVj)±{r) and the analytic formula Eq. pOf) 



quite well, except for a Ict discrepancy at separations below 30 Mpc/h. The error range is the same as for Figure [2] 
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FIG. 4: The solid red line shows vfj^ computed from the Virgo simulation using only the radial velocities, Eq. (|15p . while the 
dashed line shows Vij and shaded la errors computed using all three velocity components, the same as in Fig. [2] 



Measurement of the radial velocity of individual clusters via their kinematic Sunyaev-Zeldovich signal is affected by 
various error sources, including detector noise in the microwave maps, separating the small signal from other larger 
signals at the same frequencies (particularly the thermal SZ signal, infrared point sources, and gravitational lensing by 
the cluster), the internal velocity dispersion of the intracluster medium, and X-ray temperature measurement errors. 
In this Section, we call the total error from all of these sources "measurement error." We also consider separately the 
errors arising from cosmic variance and Poisson noise; both of these error sources are independent of the measurement 
errors for any individual cluster. 



Upcoming multi-frequency Sunyacv-Zeldovich measurements with arcminute resolution and few /iK sensitivity 
have the potential to obtain galaxy cluster peculiar velocities. However, the kinematic Sunyaev-Zeldovich signal is 
small compared to the thermal SZ signal, and is spectrally indistinguishable from the primary microwave blackbody 
fluctuations or their gravitational lensing. In addition, radio and infrared galaxies contribute substantial signal in the 
microwave bands, and are expected to be spatially correlated with galaxy cluster positions [i^. Comparatively modest 
error sources can substantially hinder cluster velocity measurements if they are not well understood and accounted 



Major potential sources of error in measuring the velocities of individual galaxy clusters include internal cluster gas 
velocities, the confusion-limited noise from point sources, uncertainties in extrapolating measured point sources to 
the frequencies of a particular experiment, instrumental noise, and the particular frequency bands available. Previous 
studies shows that primary microwave background fluctuations plus point sources set a confusion limited velocity 
error of around 200 km/s for an experiment with arcminute resolution and few /iK sensitivity [24, 49, 50], provided no 
other point source follow-up observations are utilized. The bulk flow of the gas in the intracluster medium contributes 
to an irreducible error of 100 to 150 km/s Also, Ref. [l^l shows that to extract velocity from SZ observations 

at the three ACT frequency channels (145, 220, and 280 GHz), a followup measurement of X-ray temperature of 
the cluster is needed to break a spectrum degeneracy between cluster gas velocity, optical depth, and temperature. 
While Ref. [2^ studied over 100 simulated clusters, the rest of these studies use only a few. All of these error sources 
require detailed simulations of particular experiments observing realistic simulated clusters and optimal algorithms for 



IV. ERROR SOURCES 



A. Velocity Measurement Errors 



for. 
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extracting cluster velocities from measurements in particular frequency bands and at given instrumental noise levels. 
The ultimate distribution of velocity errors is still uncertain and future study in this direction is needed. In order to 
study the effect of measurement errors on parameter estimation, we make the simple assumption that velocity errors 
have a normal distribution with a magnitude between 100 and 500 km/s. Directly adding all of the known sources 
of error from previous studies gives velocity measurement errors typically in the range of 400 to 500 km/s; however, 
with further imdcrstanding of systematic errors and point sources, the error budget may be reduced. 



In addition to cluster velocity, we must measure cluster redshift to construct the estimators of the mean pairwise 
velocity and the velocity correlation, which involve knowledge of the separation vector between the two clusters. 
For clusters at cosmological distances, the Hubble contribution to its redshift will typically be much larger than its 
peculiar velocity contribution, which we can also correct for with a direct velocity measurement, so direct error in 
the cluster redshift will be the largest contributor to the cluster position error. Typically, we will be concerned with 
cluster separations larger than 30 Mpc/h, for which the cluster velocity field is in the mildly nonlinear regime and 
can be well described by the halo model approximation. 

A redshift error of 500 km/sec corresponds to a direct Hubble distance error of around 5 Mpc/h, typically only 25% 
of the closest cluster separation of interest; even for redshift errors of 1000 km/sec, most pair separations will not 
be dominated by this error. For the remainder of this paper, we assume that the cluster sample for which velocities 
are determined also have spectroscopic redshifts from which their distances are determined, and we assume that the 
distance error effect on the cosmological parameters will be negligible compared to the direct velocity errors. For 
spectroscopic measurements of many galaxy clusters, the distance to lowest order is simply determined by the average 
of the galaxy redshifts, with an error given roughly by the cluster galaxy velocity dispersion divided by the square 
root of the number of clusters' galaxies. Cluster line-of-sight velocity dispersions will typically be 500 km/sec, so 
multi-object spectroscopy can clearly provide adequate redshift measurements. The systematic error induced because 
not all clusters will be virialized is potentially important, although beyond the scope of this paper. 

Spectroscopic redshifts for a galaxy cluster at z = 1 requires roughly an hour of observation on an 8-m class telescope. 
Spectroscopic follow-up of hundreds of clusters per year is a large program for a single telescope; spectroscopic redshifts 
for thousands of clusters will comprise a multi-year program on more than one telescope. This is likely to be a 
significant portion of the effort and expense in building a cluster peculiar velocity survey with thousands of clusters. 
Note that cluster galaxy spectroscopic redshifts are also valuable for dynamical mass estimates; see, e.g., [5l|, [s^ . 
The ACT collaboration has plans for spectroscopic follow-up observations of SZ-detected clusters using the Southern 
African Large Telescope (SALT), a new 10- meter class instrument. If only photometric redshifts are available, typically 
giving a distance accuracy of one to two percent times 1 + z, cosmological constraints must be re-evaluated. In general, 
constraints will be less stringent, although it is not immediately clear whether the resulting distance errors will have 
an effect which is significant compared to the velocity errors. In our case, redshift errors propagate only into the 
geometric portions of the mean pairwise streaming velocity and velocity correlation estimators, but the velocity errors 
are unaffected. This issue will be addressed in detail elsewhere. 



In addition to measurement errors for individual cluster velocities, cosmological quantities are also subject to errors 
from cosmic variance (any particular region observed may have different statistical properties from the average of the 
entire universe) and Poisson errors due to the finite size of the cluster velocity sample used to estimate the velocity 
statistics. Here we discuss these errors for each of the three velocity statistics. Detailed derivations of the expressions 
in the rest of this Section are given in the Appendices. 



Consider a cluster velocity survey with a measured redshift for each cluster. For the probability density function, we 
write cosmic covariance between two different velocity-redshift bins [v, z]i and [v, z]j as C^"^, which can be expressed 
as 



B. Redshift Errors 



C. 



Cosmic Variance and Poisson Noise 



1. Probability Density Function 
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FIG. 5: The effect of measurement errors on the velocity probabihty distribution function: from top to bottom, velocity 
measurement errors of cjv =100, 200, 300, 500, and 1000 km/s. Also shown are the probability distribution function evaluated 
directly using the Virgo lightcone numerical simulation (dotted curve with error bars) from Figure [T] 



where 

f d'mmb('m,a)n(m)p(v\m,S,a) 

nv[v, z) = — (24) 

J avarnn[m) 

and is the comoving length of the redshift bin within the sky survey region [53| . 

For Poisson errors, let be the total number of clusters in bin i. We are interested in the error in rii — Ni/N^ with 
Nz the total number of clusters in a particular redshift bin summed over all velocities; the measured corresponds to 
the theoretical quantity ny{v, z), Eq. PH)) . integrated over the velocity-redshift bin [v, z\i. The expression for Poisson 
errors can be written as 



5n, = (V^, + n,)/^/Nz (25) 

where the first term is from the error in Ni and the second from the error in N^- 

Random velocity measurement errors will smear out the velocity PDF. We quantify the effect of measurement errors 
by convolving the PDF with a normal distribution of velocity errors. 



n 



obs / 



{v,6v,z)= dv dv'f{v\z)exp[~{v'~v)^/2al] (26) 



where tr^ is the dispersion of the normally distributed velocity errors and the integral is over the velocity bin. Then 
the expression for the total covariance can be written as 

CT" {v,,zf, v,,z^) - C"" (zj) + (<5n,)25,, (27) 

The various curves in Fig. [5] show the effect of random velocity errors of different sizes, Eq. (j26|l . while the top 
dotted curve with shaded error region gives the actual value for the probability distribution function from the VIRGO 
simulation with Poisson plus cosmic variance errors. Smearing the distribution by random velocity errors is largely 
degenerate with the effect of varying cosmological parameters. This means that the velocity probability distribution 
function as a probe of cosmology is limited by how well the measurement error can be understood from simulated 
measurements. 
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FIG. 6: Fractional errors 5vij /vij for a cluster velocity survey covering 5000 square degrees: the red square points represents 
the Poisson error; black triangles represents cosmic variance and the Blue lines represents measurement errors (from bottom 
to top ai, = 100, 200, 300, 500 and 1000 km/s). Note that all the errors scales as yfj aky for other survey areas. 



2. Mean Pairwise Streaming Velocity 

The mean pairwise streaming velocity statistic is binned in pair separation and redshift. The cosmic covariance 
between two bins [r, z]p and [r, z]q can be written as 



We add in quadrature the Poisson error and measurement error for ripair cluster pairs and write the total covariance 

as 

C-^ (r„ r„ z,) = Cl^^M + { + ^ ) ^p'^ (29) 

\ '^pair '^pair / 

Figure [6] plots fractional errors for Vij as a function of pair separation for a survey area of 5000 deg^. For a survey 
area /sky, fractional errors scales as roughly -^//sky. Note that the Poisson error decreases for larger separation since 
more clusters pairs are available to average over, whereas cosmic variance has an increasing effect at larger separation. 
The combined effect of cosmic variance plus Poisson errors dominates the error budget when velocity measurement 
errors are below 200 km/s. Note that even when the measurement errors are as high as = 500 km/s, the total 
error is typically 50% of the magnitude of mean pairwise streaming velocity. We will show in Sec. |V] that this fact 
makes mean pairwise streaming velocity a potentially useful probe to study cosmology. 



3. Velocity Correlation Function 



Similarly for the velocity correlation function, the expression for cosmic covariance can be written as 



C. 



(ViVj) 



Stt 



dlnDr, 



din a 



din Da 
din a 



X / dk3i{krp)j^{kr,)[P(k)Y{p)i,{q); 



(30) 
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FIG. 7: Same as in Figure |6]for the fractional error 5{l^ViVj)) / (viVj) . 



using the notational abbreviation 



J dm (T_i 



(Jo (to) 



-i{m) 



(31) 



In Eq. ([50)1 we have ignored the contribution of the second {I2) term in Eq. (|20p . At larger separations relevant 
here, this term, being weighted by (,{r), is an order of magnitude smaller than the first term and hence has negligible 
contribution to the cosmic variance. 

Again we add in quadrature the Poisson error and measurement error for ripair cluster pairs and write the total 
covariance as 



{v,v.j){r,z) 



1 



(32) 



Figure [7] shows the various errors in the velocity correlation function. The trends are similar to those for mean 
pairwise streaming velocity. Measurement errors dominate the error budget for > 200 km/s. Note however 
the increase in fractional errors with the increase in measurement errors. For ay = 500 km/s, the contribution of 
measurement errors to the total error is almost 90%, nearly double that for the case of mean pairwise streaming 
velocity. 



V. CONSTRAINTS ON DARK ENERGY PARAMETERS 



Now we consider constraints on dark energy parameters for various survey areas and over a range of velocity errors. 
Following the Dark Energy Task Force, we describe the dark energy in terms of three phenomenological parameters: 
its current energy density JIa, and two parameters wq and Wa describing the redshift evolution of its equation of state 
w{a) = Wq + {1 — a)wa. Assuming a spatially flat universe, the set of cosmological parameters p on which the velocity 
field depends are the normalization of the matter power spectrum as (or equivalently the normalization constant B in 
Eq. ([5])), the power law index of the primordial power spectrum ng, and the Hubble parameter /i, plus the dark energy 
parameters. We perform a simple Fisher matrix analysis to find constraints on these parameters from measurements 
of the three velocity statistics described in Sec. HTl 

We consider a fiducial model similar to that assumed in the DETF report [l3| with erg = 0.9, = 1, /i = 0.7, 
JIa = 0.72, Wq = —1, Wa = 0. To make quantitative comparisons with the conclusions of the DETF report, we 
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FIG. 8: The change in Ict parameter constraints with velocity error (normal distribution of width (t„) for a 4000 deg^ survey 
area, for the three statistics n„ (blue dashed), Vij (red short dashed) and {viVj) (black solid). The four panels are for the 
parameters (top left), SIa (top right), Wa (bottom left), and the Figure of Merit (bottom right). 



compute values for the expression [o'(wo)cr(u'p)]~^, which is listed in the DETF summary tables. We refer to this as 
the "Figure of Merit" (FOM) for convenience, although this term refers to a slightly different quantity (inverse area 
of the ellipse of 95% confidence limit in the Wp — Wa plane) in the DETF report. Here Wp is the equation of state at 
the pivot point defined as Wp = wq + {1 — ap)wa with Up = 1 + [F^^]wQWa/[F^^]waWa ^^'^ P the Fisher information 
matrix for a given experiment. 
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TABLE I: la errors on dark energy parameters for a 4000 deg^ survey area plus cosmological priors from Planck and HST 
0, E^li assuming a spatially flat cosmology. 
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(ViVj) 


Vij 
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100 


0.06 


0.083 


0.099 


0.16 


0.26 


0.2 


0.007 


0.007 0.016 
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94 


200 


0.1 


0.1 


0.11 


0.29 


0.33 


0.25 


0.014 


0.008 0.018 


60 


76 


71 


300 


0.18 
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0.14 


0.53 


0.42 


0.34 


0.026 


0.009 0.019 


20 


54 
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500 


0.39 


0.18 


0.25 


1.32 


0.61 


0.65 


0.046 


0.012 0.026 


5 


31.5 


21 


1000 


1.28 


0.31 


0.9 


4.7 


1.11 


3.0 


0.060 


0.018 0.048 
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14.5 


3.0 



TABLE IL Same as Table HH for a 2000 deg^ survey area. 
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0.25 
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1.48 


4.2 
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0.022 
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1.6 



The Fisher information matrix for each of the three statistics is [SJl 

-Q^ (33) 

where stands for either n^, Vij{r, z) or {viVj){r, z), C'^{ij) is the total covariance matrix in each bin for the statistic 
(j), Eqs. (P7)) . (|29p . and ([5^ . and the partial derivatives are evaluated for the fiducial values of the cosmological 
parameters. The values i and j index the bins [r, z]i and [r, z]j for the mean pairwise streaming velocity and velocity 
correlation function, while for (j> — iiy, i and j refer to [v, z]i and [v, z]j. The inverse of the Fisher matrix has diagonal 
elements which are estimates for the variances of each cosmological parameter marginalized over the values of the 
other parameters, and the non-diagonal elements give the correlations between parameters. 

Figure [8] shows the degradation of parameter constraints with increasing velocity error ay for a 4000 deg^ survey 
area. It is evident that parameter constraints from Vij are more robust to increases in velocity error than those from 
n„ and (viVj). This is because Suij depends linearly on cr^, while 6{viVj) varies as and for the distribution gets 
smeared with increases in CTu. Constraints on wq, Wa and change roughly by a factor of two and the constraint 
on the FOM by a factor of three, for the factor of five increase in ay from 200 to 500 km/s. Compare this to the 
corresponding change for {viVj): wq, Wa and flA constraints change roughly by a factor of 6 to 8 and the FOM 
constraint by a factor of 30 for a similar change in a^. For n^,, the corresponding degradation in constraints are 



TABLE III: Same as Table HI for a 400 deg^ survey area. 







Wo 






Wa 












[(j(wo)o-(wp)] ^ 






(ViVj) 


Vij 


n„ 


(ViVj) 


Vrj 


n„ 


{ViVj) 


Vij 


n„ 


{viVj) 


Vij 


n„ 


100 


0.13 


0.20 


0.2 


0.45 


0.72 


0.51 


0.019 


0.015 


0.023 


30 


22 


29 


200 


0.24 


0.25 


0.24 


0.76 


0.92 


0.64 


0.034 


0.017 


0.026 


11 


16 


21 


300 


0.41 


0.31 


0.31 


1.39 


1.15 


0.85 


0.048 


0.020 


0.031 


4.0 


11 


14 


500 


0.92 


0.6 


0.53 


3.4 


1.66 


1.53 


0.058 


0.024 


0.044 


1.4 


0.7 


5.2 


1000 


3.6 


0.78 


2.42 


13.3 


3.0 


8.0 


0.061 


0.033 


0.061 


0.38 


3.3 


0.7 



15 



TABLE IV: 1(7 constraints for dark energy parameters for a fiducial cluster survey of 4000 clusters with velocity errors cr„ = 1000 
km/s, for cluster number counts, cluster velocities, and the two combined. Planck and HST cosmological priors [13: [13 Pl^^ 
spatially flat cosmology assumed. 



Parameters Priors Counts Velocity Combined 



nA[o.7] 


0.062 0.052 


0.033 


0.025 


wol-l] 


0.94 


0.78 


0.52 


Wa[0] 


2.95 


3.0 


1.8 


FOM 


2.8 


3.0 


7.0 



roughly by a factor 1.5 to 3 for wq, Wa and f^A while the FOM constraint degrades by roughly a factor of 4. Table |T] 
lists the constraints as a function of velocity error for a 4000 deg^ survey area, while Tables [TTl and lllll give constraints 
for 2000 deg^ and 400 deg^ respectively. 

Note that the velocity correlation function (viVj) provides the best constraints on the dark energy equation of 
state {wq, Wa, and FOM) for (T„ < 200 km/s. It might be possible to achieve such values of errors in future surveys 
with better understanding of point source contamination and other systematics. However for more realistic near-term 
errors of 500 km/s, the mean pairwise peculiar velocity provides better constraints on dark energy parameters, 
and this statistic will be used in the following sections which consider how cosmological constraints will be improved 
by using cluster velocity information. 

VI. COMPLEMENTARITY OF CLUSTER VELOCITIES WITH CLUSTER NUMBER COUNTS 

For a given SZ survey, we can potentially obtain both cluster counts and cluster peculiar velocities. Given these 
two different data sources from the same survey, what is the joint constraint on dark energy parameters they provide? 
Consider a fiducial Stage II survey of 4000 galaxy clusters proposed by the DETF report [l3| (see Table |V] for 
details), plus the addition of cluster velocities with measurement error ay = 1000 km/s, along with cosmic variance 
and Poisson errors to estimate the mean pairwise peculiar velocity statistic Vij. This is not a particularly stringent 
velocity error, and it is likely obtainable with currently planned surveys with foreseeable follow-up observations or 
theoretical assumptions about cluster properties. Table HVl gives the constraint on the dark energy parameters derived 
considering cluster counts only, considering cluster velocities only, and the joint constraint from both. Also given are 
HST plus Planck prior constraints assuming a flat spatial geometry. We find cluster velocities provide a better 
constraint on f^A and wq than cluster counts, even for a measurement error of ay = 1000 km/s. The constraint on 
Wa is comparable for the two probes. The combined constraint is a factor of two better than the counts-only case for 
Oa, Wq and the Figure of Merit, and at least a 60% improvement for Wa- The relative complementarity between the 
two probes is shown in Figure [5) 

We have assumed that the cluster velocity and cluster density observables are statistically uncorrelated. As they will 
likely be obtained from the same set of clusters, it is reasonable to ask whether this is actually true. A straightforward 
analytic calculation shows that the cross-correlation between velocity and density will be proportional to the matter 
bispectrum, so we expect it to be small compared to the signal from the velocity correlations, which are proportional 
to the matter power spectrum. We intend to confirm this prediction from sets of large-volume numerical simulations 
when these are available. 



VII. COMPARISON WITH DETF PROPOSED EXPERIMENTS 

The Dark Energy Task Force report [l3| considers four different potential probes to study dark energy parameters: 
weak lensing(WL), baryon acoustic oscillations (BAO), cluster counts (CL) and SNIa (SN) luminosity distance mea- 
surements. The relative merits of these probes have been discussed in detail in the DETF report both for ongoing 
and future projects. In this section we compare our fiducial velocity survey with each of the four DETF probes. To 
assess the advantage of adding cluster mean pairwise peculiar velocity Vij as a dark energy probe, we have considered 
only the most optimistic forecasts for the DETF surveys (i.e. survey assumptions that provide maximum constraint 
to the FOM assuming a flat universe plus HST and Planck priors) for each Stage in the DETF report. Table fVl gives 
a brief description of the DETF surveys considered here and our corresponding assumed cluster velocity surveys. We 
have used the actual Fisher matrices used by the DETF team along with their priors for the following comparisons. 
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FIG. 9: The relative complementarity of velocity and cluster counts. Shown are la error ellipses in the wo — Ha plane (left) 
and the Wa — plane (right) for 4000 clusters with normally-distributed velocity errors of ~ lOOOfcm/s. The three ellipses 
are for cluster velocities (red), cluster counts (blue) and the combination of both (black). Planck and HST cosmological priors 
[3, \EM ^^'^ ^ spatially flat cosmology are assumed. 



TABLE V: Parameters defining various surveys discussed in the DETF report [3] , plus various cluster velocity surveys discussed 
here. 



Stages 


VEL 


WL 


SNIa 


CI BAG 


II 


N,i = 4000, /,ky = 0.01 /,ky = 0.0042 
M„,„ > 2 X lO"M0//i 
z=0. 1-1.4 


SNLS 
700 SNIa 
z=0. 1-1.0 


Nci = 4000 None 
/sky = 0.005 


III 


Nci = 15000 

/sky = 0.05 


DES 

/sky = 0.1 


2000 SNIa 
Spectroscopy 


Nci = 30000 /sky =0.1 


IV 


Nci = 30000 

/sky = 0.1 


SKA-o 

/sky =0.5 

z = 0.1-1.7 


Space 
2000 SNIa 


Nci = 30000 SKA-o 

/sky = 0.5 /sky = 0.5 

z = 0-1.5 



For the fiducial cluster velocity surveys, we have assumed that SZ surveys will be sensitive enough to detect the 
kSZ signal from all clusters with M > 2 x IO^'^Mq/H. To be consistent with the DETF report, the total number of 
clusters for each survey corresponds to as ~ 0.9. If cts — 0.76 |2i| is used, then the corresponding number of clusters 
decreases by a factor of 30%. However, a velocity survey is sensitive to only the number of detected clusters and not 
the volume of the survey. So our conclusions will still be valid if the survey area is increased to compensate for a 
lower value of erg. 

A comparison of velocity with other probes is shown in Figure fTOl HST and Planck cosmological priors and spatially 
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FIG. 10: A comparison of the error in the dark energy density 6Qa and the dark energy figure of merit obtained from velocity 
statistics with that from DETF probes. The top two panels are for Stage II experiments; the dark region shows the range in 
the parameter error for the DETF- assumed ranges in the measurement errors. For cluster velocities we assume a range from 
ffv = 200 to 1000 km/sec. The middle panels show the results for Stage IV measurements. The bottom panels show the relative 
improvement in parameter measurements at Stage IV when cluster velocities are combined with all of the other DETF probes. 
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flat cosmology are assumed for afl the probes. Each plot shows a range of parameter errors for each experiment, 
corresponding to cluster velocity measurement errors ranging between 200 and 1000 km/sec, and other measurement 
errors as in the DETF report. At Stage II, velocity provides a competitive constraint on n\ compared to SNIa, and 
much better constraints than weak lensing or cluster number counts. Even a modest velocity survey would yield a 
factor of two better constrain on J^a than cluster counts or weak lensing. Cluster velocities also provide two to three 
times better constraints to the figure of merit compared to weak lensing or cluster counts at Stage II. Ultimately at 
Stage IV, however, weak lensing provides the most accurate measurements of dark energy density and the figure of 
merit. But constraints from velocity are competitive with those from supernovae and better than those from cluster 
counts or baryon acoustic oscillations. Stage II and III experiments yield an average 20% improvement in cosmological 
parameter determination, and Stage IV about a 7% improvement, when velocity information is combined with the rest 
of the dark energy experiment results. This corresponds to an improvement by factors of 1.5 to 2.5 in the dark energy 
figure of merit. These types of statistical comparisons of course assume zero systematic errors; cluster velocities will 
ultimately be more valuable than these numbers indicate, due to their completely different systematic errors from the 
other challenging techniques. All of these methods will in the end be dominated by systematic, not statistical, errors. 

VIII. DISCUSSION 

The various studies of galaxy cluster peculiar velocities in this paper yield a number of interesting conclusions. The 
measurement of peculiar velocities of objects at cosmological distances is of fundamental importance, as it directly 
probes the evolution of the gravitational potential. The kinematic Sunyaev-Zeldovich effect in clusters of galaxies 
promises a direct tracer of this signal, with errors largely independent of cluster redshift. Although the current 
uncertainty in velocity measurements is large with w 1000 km/s [55| for individual clusters, upcoming multi-band 
experiments like ACT [H^ or SPT ^] with arcminute resolution and few /iK sensitivity have the potential to measure 
peculiar velocities with velocity errors of a few hundred km/s for large samples of clusters, opening a new window on 
the evolution of the universe. We have considered three separate cluster velocity statistics here, computing them using 
the halo model and comparing with numerical results. For surveys with thousands of cluster velocities with errors 
of a few hundred km/sec, dark energy constraints competitive with other major techniques (cluster number counts, 
baryon acoustic oscillations, supernova redshift- distance measurements, and weak lensing) can be obtained from the 
mean pairwise peculiar velocity Vij, with different systematic errors. Even for velocity errors as large as 1000 km/s 
for individual clusters, a velocity catalog for several thousand clusters can improve dark energy constraints from the 
corresponding cluster number counts by a factor of two. 

Throughout this work, we have simply assumed that cluster velocities can be extracted with normal errors and 
no bias from a Sunyaev-Zeldovich sky survey of sufficient angular resolution and low enough noise. Connecting the 
measured SZ signal to the cluster velocity is a non-trivial task. The three ACT measured frequencies at 145, 220, 
and 270 GHz have a degeneracy which prevents the cluster velocity from being determined uniquely along with the 
cluster optical depth and temperature [2^.[57|. This can be remedied several ways, including adding other microwave 
bands |57i | or X-ray temperature measurements [23 | , or assuming cluster scaling relations between various measurable 
quantities [H, . 

Further complications arise because the measured signal is not due only to the Sunyaev-Zeldovich distortions, but 
also contains the blackbody primordial microwave fluctuations, gravitational lensing of the microwave background, 
infrared and radio point sources which can be correlated with galaxy cluster positions, and galactic dust (see [63| for 
a description of sky simulations incorporating all of these signals). The kinematic SZ signal must be separated from 
all of the others via a combination of frequency and spatial filtering. With sufficient data, this can clearly be done 
uniquely, but with limited wave bands, spatial resolution, and noise levels, any kSZ signal extraction will be subject 
to some amount of measurement error. Evaluation of this error for various observing strategies is important and we 
are currently pursuing it using simulations. Even with perfect separation, internal cluster gas motions provide an 
irredicible error fioor for kSZ cluster velocity measurements of around 100 km/sec [23l. [25j. 

Component separation and other issues may also lead to systematic errors. We are currently modeling systematic 
errors in velocity measurements in some detail, but it is clear that at minimum, cosmological constraints based on 
cluster velocity measurements are much less prone to systematic errors due to uncertainties in the relation between SZ 
distortion and cluster mass than constraints based on cluster number counts [isl . [6l| . This is potentially the dominant 
systematic error for cluster number count studies, and largely mitigating it is one strong incentive for pursuing cluster 
velocities as an alternative probe of dark energy. An additional advantage of using cluster velocities is that the 
cluster velocity distribution function riy should be symmetric with respect to positive and negative peculiar velocities, 
by homogeneity of the universe. Departures from symmetry are easily diagnosed and can be used as a monitor of 
unknown systematic errors. The downside of cluster velocities is that the kSZ signal is much smaller than the thermal 
SZ signal, on the order of 5 to 10 for large clusters with typical peculiar velocities. Separating this small signal 
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from other larger ones may lead to different systematic errors. But potential constraints on dark energy from cluster 
velocities are good enough, and the other methods of measuring dark energy properties are hard enough, that building 
a cluster velocity catalog with a different set of systematic errors from other techniques is surely valuable. 

A number of further lines of work related to cluster velocities are worth pursuing. Here we have considered 
three different galaxy cluster velocity statistics: the velocity probability distribution function n„, the mean pairwise 
velocity dispersion Vij, and the velocity correlation function (viVj). Each constrains well a different set of cosmological 
quantities. We have not attempted a joint analysis, finding the combined cosmological constraints from all three 
statistics: the correlations between the statistics are complicated, and no clear way to derive them analytically 
presents itself. Proper joint constraints will require numerical evaluation of the correlations between statistics from 
sets of large cosmological simulations, which is feasible but demanding. A related question is the extent to which 
these three statistics, which are convenient from a theoretical and observational point of view, exhaust the useful 
cosmological information on dark energy constraints: are there other velocity statistics which, when combined with 
these three using the correct correlations, would further tighten the constraints? This is an open, and challenging, 
question. 

On the numerical front, we have performed limited tests comparing the VIRGO simulation results with the halo- 
model expressions for the velocity statistics here, finding reasonable agreement for the particular cosmological model 
the simulation is based on. This is encouraging, but it would be reassuring to have explicit comparisons between 
theory and simulation for a wider range of models. Such computations require cosmological simulations over very large 
volumes, to capture a sufficient number of clusters with large enough masses, but can be done with fairly low mass 
resolution, since we only care about bulk cluster properties and not internal cluster details. Sets of such simulations 
are currently in progress. 

The kinematic SZ signal does not directly measure cluster peculiar velocity, but rather is proportional to a line- 

of-sight integral of the cluster gas' local peculiar velocity times its local density. Thus the kSZ effect is actually 
proportional to the cluster gas momentum with respect to the cosmic rest frame. We can sidestep the entire difficult 
observational issue of inferring cluster velocities from kSZ measurements by using cluster momenta instead. Wc then 
need theoretical calculations for the cluster momentum statistics corresponding to the velocity statistics considered 
here. Momentum statistics have the possibility of being just as cosmologically constraining, but easier to compare with 
observations. We have not found any suitable analytic approximations to the cluster momentum statistics, but this 
could also be evaluated numerically using large-volume, low- resolution N-body simulations mentioned above. The 
other related issue is connecting the cluster mass, which is used to evaluate cluster momenta in an N-body simulation, 
to the cluster gas mass, which gives the SZ signal. Wc need to understand the extent to which the cluster gas fraction 
is constant, or the extent to which we can understand its statistical distribution. We have already made initial steps 
to investigate this issue, finding, among other things, that the gas fraction in galaxy groups is affected non-negligably 
by quasar feedback, which heats the gas and suppresses star formation. However, at mass scales substantially below 
galaxy clusters, the gas fraction appears to be reasonably independent of mass. Probing this relation for clusters 
is a challenging computational issue, requiring sophisticated hydrodynamical simulations in much larger volumes to 
obtain information about galaxy clusters large enough to be of SZ interest. 

As with so many cosmological sources of information, the advent of the dark energy era has given a new urgency 
to precision measurements. Galaxy cluster velocities, obtained via their kinematic Sunyaev-Zeldovich signal, directly 
probe the growth of structure in the universe via gravitational instability. The signals are small, but the advantages 
manifest. We firmly advocate that cluster velocities should be added to the arsenal of tactics now trained on the dark 
energy issue. 
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APPENDIX A: ERRORS FOR THE PROBABILITY DENSITY FUNCTION 

1. Poisson Error 

Let Nz be the number of halos in redshift bin z + Jz, and be the number of hales in both the redshift bin z + 5z 
and the velocity bin v + 5v. In a given velocity bin, the fractional density The observable in the normalized histogram 
of cluster velocities in a given redshift bin is then = Ny/Nz- Thus ri„ suffers from uncertainties in both numerator 
and denominator. We write the uncertainty in as 

Assuming Poisson errors, SN^, = \fN\, and 5Nz = \/Nz- We write 



5n, = V^J^/Nz + nJ^/Nz - ^'"^^^^'"^ (A2) 



2. Cosmic Variance Error 

Write the cosmic covariance between two different bins [f^jZi] and [ Cj"" ; here Vi denotes a particular 

velocity bin at an epoch of redshift Zi . C^"" is defined as 

C"" (u) {{fivi - nyi){nyj - riyj)) (A3) 

where iiy denotes the estimated PDF and rim = n-u{vi,Zi) etc. Using n(r7i, (5, x) = (1 + b{'m)S{'x.))h{'m) and hy — 
V{r)^^[J cPx J dm mn{m\S,:x.)p(v\m,6)]/ J dmmn^mlS), 

C""(«j) = b{mi,Zi)b{mj,Zj)fifj (5(xi, Zi)^*(xj, z^)) (A4) 

where (...) denotes the ensemble average over the survey volume Vh and can be written as 

{5,5j)^^ j ^ d^r J J d3xd3x'M^(x)W^(x')(5(x,a)(5(x',a')<5|5(x-x'-r) (A5) 

where 

5(x, a) = DaS{x) j d3k(5(k)e*'^ '', (A6) 

W{x) is the tophat window function defined after Eq. ([T]) and (5(x) is the field describing linear comoving density 
perturbations evolved to the present; the three-dimensional Dirac delta distribution is written as (5|j(x). We can then 
write 

{5^5^) = !^^£^ f ^3j. f f d^^d^x'W{x)W{x')S{ji)S{^')Sl{^-^' -r) 
Vn J J 

where we write conventionally [6^ 



J J J d^kd^k'5{k)S*ik')e-'^"h{k~k',r). (A7) 



h(k, r) = / d^xW(x)W(\x + rDe'*"^. (A8) 

V[r) J 

In the limit of a survey region large compared to the scale r, h{k,r) ^ S'^{k), r ^ Rq [g^, [6^ with the convenient 
notation Vh = ^ttRq/S for a spherical survey volume, giving 

J d^xWix)Wi\x + r\)e'^^~^'^-'' J d^xW^x) cx ^^^(k - k'). (A9) 
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Then 



so C"" {ij) can be written as 



where 



Vh 



J dkPP{k)jiikRn), 



Rn 



-n,nj J dkk^P{k)jx(kRn) 



z) 



J dmmb{m,a)n(m)p(v\m,S,a) 
J dmmn{'m) 

which is equivalent to Eq. (j23p . The expression p(w|r7i, (5) is defined in Eq. ([S]). 



(AlO) 



(All) 



(A12) 



APPENDIX B: ERRORS FOR THE MEAN PAIRWISE STREAMING VELOCITY 

1. Poisson Error and Measurement Error 

We begin with Eq. (|15p for the estimator of the mean pairwise streaming velocity. Assume a particular velocity is 
measured with an accuracy Sv. So the error 6vij in Vij can be written as 



so that 



— f V2(j.„ + Vi 



(Bl) 



(B2) 



where we have used drip ~ \fnp assuming a Poisson distribution, and 

5Sy K - wj] ^ \/2[(5t;2 + + ... + ji/2 ^ V2^^pa„. (B3) 
Here the individual velocity errors are added in quadrature and the last line follows from the central limit theorem. 



2. Cosmic Variance Error 



The cosmic covariance for mean pairwise streaming velocity between two separation and redshift bins \rp , Zp] and 
\rq , Zg] can be written as 



C"''^ (pg) = {(v.,j (p) - v,j (p)) {v,j (q) ~ v,j (q))> = (wy {p)vij (q)) ~ v^j {p)v.^ (q) 



(B4) 



where Vij is the estimated mean pairwise streaming velocity from the survey volume and Vij is its cosmic mean value, 
(vij) = Vij. Using the expression for mean pairwise streaming velocity given in Eq. (jlip . the above expression can be 
written as 



_rpH{ap)ap{^^ 



2 „, . f dlnPg 



where ^ is an estimator for the volume-averaged correlation function 

1 



V{r) Jo 



dr' r'2^(r') 



(B5) 



(B6) 
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An estimator ^(r) for the two-point correlation function ^(r) is 



|(r) = -1- / d^x'Wix') I d^xWix)5{x)6ix')6l{^ - x' - r), (B7) 



V{r 

so an estimator for the volume-averaged correlation function can be written as 

If , 1 

'v(r 



ar) -vV)Jvr '^'^'vW) J ^'^^(^) J d'^'Wix')S{m^')Sli^ - x' - r') (B8) 



where the survey volume is given by ^(r) = Jd'^xVl^(a;)T4^(|x-|-r|) for a normalized window function J (PxW{x) = 1. 
Fourier transforming S{x), we can write 

|(r) = -1^ / d^r'-^ f dh' f d'^jiWix) f d'^^'W{x')Sl{^ - ^' - r') [ [ d3kd=^k'(5(k)r (k')e'('' "-''' "') 

Jv(r) y{f') J J J J J 

= I dh' f f dVk'6{k)S*Ck')e-'^ '''h(k^k\r'). (B9) 

V(r) Jo J J 

Using {^{r)) = ^(r), we can then write 

= , , , , / d^re-'^-^hCk-k^r) / d^r' e''^ " h* (k - k' ,r') 
V{rp)V{rq) Jo Jq 

X J d'^kj d'^k' J d^k^ j d^k; [(5(k)(5*(k')<5(ki)<5*(k'i)) - ((5(k)r(k')) (5(ki)r (k'^))] . (BIO) 

The term in brackets can be written as 

[...] = 6l{k + k^)P{k)5l{k! + k'i)F(fc') + (533(k - k'i)F(fc)5|5(k' - ki)F(fc'). (Bll) 
Substituting this expression into Eq. (|B10p gives 

&^{pq) = , } , , I rf^k / d3k'P(fc)P(fc') fe*'"'^'"'"') + e"'*' '"-^*'' '"'') / ''d3r/i(k-k',r) / " d^v'h*(k~k' ,y'). 
y{Tp)V{rq) J J V /Jo Jo 

(B12) 

As in the previous appendix, for large surveys such that r « Rn = {SYn/iTrY^^ , h{k - k',r) - Sl,{k - k') and 

_ / d^KWHx)W{\^ + r\)Wi\K + r'\) J_ 

Virp)Vir,) ^ Vn ^ ' 

So Eq. (jB12p can be written as 

cHp,) ^ ^^^n^J fj £ ^^^^^^^ 

= / dke\P(k)\^3i{krp)]i{kr,) (B14) 

ynrpTq J 

Substituting the above result in Eq (jBSp . we obtain the final expression for cosmic covariance as 



9Vn l-|-e'''°(rp,ap) 1 -he'^'°(rg,a,) V dhia J \ d\na J ' 
On scales of interest, < 1, so Eq. (|BT51) reduces to Eq. ([28]) . 
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APPENDIX C: ERRORS FOR THE VELOCITY CORRELATION FUNCTION 
1. Poisson Error and Measurement Error 



The expression for the perpendicular velocity correlation {viVj)±{r) for a particular separation r can be written as 



(CI) 



where we abbreviate = ([fi — rj] x v^) • ([ri — rj] x vj) the product of the velocity components perpendicular 

to the direction connecting the two positions. As before, is the velocity of halo i, which is measured with a normal 
error in its magnitude of dv, and Up is the number of pairs in the survey volume for a given separation distance r. For 
the rest of the appendix, we drop the perpendicular subscript for convenience. So the measurement error in (uiVj) 
can be written as 



{viVj) + S{viVj) ~ — T,ij[viV.j + 2vjSvi + 6vi5vj] 

Tip 

S{viVj) — — T,ij[2vj5vi + SviSvj] 

Tip 

= J_j:[s{v^) + {Svf] 



(C2) 



Similarly, the Poisson error is {viVj)[Snp/np] = {viVj) / ^/n^. 



Cosmic Variance Error 



The cosmic covariance for the velocity correlation function between two bins [r^, Zp ] and [^5,2^], one of separation 
Tp at epoch Zp and the other of separation Vq at rcdshift Zq, can be written as 



(C3) 



As in the case of Vij{r), we first derive an estimator for ViVj{r). In linear theory, v{k) = 5{k)/k, so v{x) 
J (Ph[d{k)/k] exp(ik • x). Then an estimator ViVj{r) measured at a separation r is 



uivjir) 



^ / d^x'W(x')vx' [ d^xW(x)v(x')Sl(x-x' -r) 
(r) J J 

/i(k-k'). 



kk' 



(C4) 



The only difference between Eq. (|C4[) and Eq. (|B9p is the added factor of kk' in the denominator. 

The expression for the velocity correlation function given in Eq. (j20p consists of two terms, expressions for which 
are given in Eqs. (j2ip and (j22p . For simplicity, here we derive the cosmic covariance of the first term using the linear 
theory expression for the velocity correlation, Eq. (jC4|l : the derivation can be easily extended to the halo model 
expression for {viVj) given in Eq. (pO)) . As argued before, the second term in Eq. (|20|) can be neglected compared to 
the first term because ^ (r) is negligible at separations of interest for r > 30 Mpc. The linear theory counterpart for 
Eq. ([^0]) can be written as 



d\nDa 

H(a)— aDa 

a In a 



SV{r) Jo 



d^kd-'k' 



,3, ,<5(k)J*(k')„_,k., 



kk' 



/i(k-k',r'). 



(C5) 



Note that this integrand is similar to that in to Eq. ([20|) . apart from the halo number density and bias factors. The 
factor of 1/3 in Eq. (jCSp . compared to Eq. (jB9|) . is because only the radial velocity components are considered. 
Proceeding analogously to Eqs. (jB9p to (|B12p . we obtain 



C^^ipq) = ayqDl^DlH^{ap)H\aq) 



din Da 



din a 



din Da 



din a 



^ dkPikf''^^'^^''^^''^'^ 



1/2 



kVrt 



kr„ 



(C6) 



24 



This is the cosmic covariance for the hnear theory counterpart of Eq. (|20|) . Including the extra halo model factors 
gives Eq. ([5(1)) . 
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